Fermi-surface calculation of the anomalous Hall conductivity 
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While the intrinsic anomalous Hall conductivity is normally written in terms of an integral of the 
electronic Berry curvature over the occupied portions of the Brillouin zone, Haldane has recently 
pointed out that this quantity (or more precisely, its "non-quantized part") may alternatively be 
expressed as a Fermi-surface property. Here we present an ab-mitio approach for computing the 
anomalous Hall conductivity that takes advantage of this observation by converting the integral over 
the Fermi sea into a more efficient integral on the Fermi surface only. First, a conventional electronic- 
structure calculation is performed with spin-orbit interaction included. Maximally-localized Wannier 
functions are then constructed by a post-processing step in order to convert the ab-initio electronic 
structure around the Fermi level into a tight-binding-like form. Working in the Wannier represen- 
tation, the Brillouin zone is sampled on a large number of equally spaced parallel slices oriented 
normal to the total magnetization. On each slice, we find the intersections of the Fermi-surface 
sheets with the slice by standard contour methods, organize these into a set of closed loops, and 
compute the Berry phases of the Bloch states as they are transported around these loops. The 
anomalous Hall conductivity is proportional to the sum of the Berry phases of all the loops on all 
the slices. Illustrative calculations are performed for Fe, Co and Ni. 

PACS numbers: 71.15.Dx, 71.70.Ej, 71.18.-(-y, 71.20.Be, 75.47.-m. 



I. INTRODUCTION 

It is by now well established that the intrinsic Karplus- 
Luttinger mechanism^ plays a significant role in the 
anomalous Hall conductivity (AHC) of fcrromagnets. 
This contribution can be expressed as an integral of the 
/c-space Berry curvature over the occupied portions of 
the Brillouin zone (BZ).^"^ First-principles calculations 
of the intrinsic AHC have been carried out by several au- 
thors, using either a Kubo linear-response formula^'^ or 
a direct "geometric" evaluation of the Berry curvature,* 
and achieving good agreement with experimental values 
for several fcrromagnets. These studies revealed that the 
Berry curvature is very sharply peaked in certain regions 
of the BZ where spin-orbit splitting occurs near the Fermi 
level. As a result the calculations tend to be rather de- 
manding; in the case of bcc Fe, for example, millions of 
fc-points must be sampled to achieve convergence.^ More 
efficient approaches are therefore highly desirable. 

In a preceding paper,* we developed a strategy for cal- 
culating the AHC in which Wannier interpolation of the 
Bloch functions was used to circumvent the need to per- 
form a full first-principles calculation for every fc-point. 
Thus, while the required number of fc-points was not re- 
duced, the computational load per fc-point was greatly re- 
duced. In this approach, the actual first-principles calcu- 
lations are performed on a comparatively coarse fc-mesh. 
Then, in a postprocessing step, the calculated electronic 
structure is mapped onto an "exact tight-binding model" 
based on maximally-localized Wannier functions.^ Work- 
ing in the Wannier representation, the Berry curvature 
can then be evaluated very inexpensively at each of the 
fc-points of the fine mesh needed for accurate evaluation 
of the AHC. 



Recently, Haldane has shown that while the intrin- 
sic AHC is usually regarded as a Fermi-sea property 
of all the occupied states, it can alternatively, and in 
some ways more naturally, be regarded as a Fermi-surface 
property. ^° (More precisely, Haldane showed that these 
quantities are equal modulo the quantum of transverse 
conductivity that is well-known from the quantum Hall 
effect, since one cannot rule out the possibility that, e.g., 
some occupied bands carry non-zero Chern numbers.^") 
By a kind of integration by parts, Haldane showed how 
the integral of the Berry curvature over the occupied por- 
tions of the BZ could be manipulated first into a Fermi- 
surface integral of a Berry connection, and then ulti- 
mately into a Fermi-surface integral of a Fermi-vector- 
weighted Berry curvature, augmented with some Berry- 
phase corrections for the case of non-simply-connected 
Fermi sheets. In discussing his Eq. (23), Haldane men- 
tions in passing that this expression can also be refor- 
mulated in terms of the Berry phases of electron orbits 
circulating on the Fermi surface. 

In this paper we present a tractable and efficient com- 
putational scheme that takes the latter point of view as 
its organizing principle. In our approach, the BZ is di- 
vided into a fine mesh of equally-spaced slices normal to 
the direction of the magnetization, and the integral of the 
Berry curvature over the occupied states of a given slice 
is transformed into a sum of Berry phases of Fermi loops 
lying in that slice. As a result, the three-dimensional BZ 
integration is avoided, and the method relies instead only 
on information calculated on the two-dimensional Fermi 
surface. As in Ref. 8, an important ingredient of our ap- 
proach is the use of a Wannier interpolation scheme to 
lower the cost further by eliminating the need for a full 
first-principles evaluation at each point on the Fermi sur- 
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face. Combining these two complementary strategies, we 
arrive at a robust and efficient method for computing of 
the AHC in ferromagnetic metals. 

The paper is organized as follows. In Sec. II we 
present the necessary formulas relating Berry phases on 
the Fermi surface to the AHC, as well as their evalua- 
tion in the Wannier representation. The details of the 
first-principles calculations and the determination of the 
Fermi loops arc given in Sec. III. In Sec. IV the method 
is applied to the transition metals Fe, Co and Ni. A dis- 
cussion of issues of computational efficiency is given in 
Sec. V, followed by a brief conclusion in Sec. VI. 



II. METHOD 

A. Fermi-loop formula 

Our starting point is the AHC expressed as an antisym- 
metric Cartesian tensor in terms of the Berry curvature, 



(1) 



where the integration is over the three-dimensional BZ 
and the occupation function /n(k) restricts the sum 
to the occupied states (we work at zero temperature). 
f^n,Q/3(k) is the Berry-curvature matrix of band n, de- 
fined as 
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where M„k is the periodic part of the Bloch function V'nk- 
Because fin, a/3 is antisymmetric, we can represent it in- 
stead in axial-vector notation as 



9 ^ ^ ^a(3^^n,af3 



(3) 



a/3 



or equivalently, 0„^c«/3 = iaf3-f^n-y , where Ca/s^y is the 
antisymmetric tensor. The Berry curvature can also be 
written as 



n„(k) = Vk X A„(k) 
where the Berry connection is 

A„(k) = i(w„k|Vk|w„k) • 
Following Ref. 10, we rewrite Eq. (1) as 



O'a/3 = 



h (27r) 



X2 ^afijKnj 



(4) 
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where 



^[ dk/„(k)n„(k). (7) 



For the case of a completely filled band lying entirely 
below the Fermi level, Haldanc has shown^" that K„ is 
quantized to be a reciprocal lattice vector (the "Chern 
vector" ) , as will become clear in Sec. II B below. 

Let sii and bj be a conjugate set of primitive real-space 
and reciprocal-space lattice vectors respectively, • bj = 
2-K5ij, and let 



so that 



K„ — Cnj bj 



(8) 



(9) 



In order to compute c„3, for example, we choose the BZ 
to be a prism whose base is spanned by bi and b2 and 
whose height is 27r/a3, and convert the integral into one 
over slices parallel to the base. In general, separate cal- 
culations in which the slices are constructed parallel to 
the b2-b3 and bi-ba planes are needed to compute c„i 
and c„2 respectively.-'^-^ However, this can be avoided in 
the common case that the magnetization lies parallel to 
a symmetry axis; one can then choose bi and b2 perpen- 
dicular to this axis, and only c„3 needs to be computed. 
Inserting Eq. (7) into Eq. (8) yields 



(t>n{kA_) 



where 



'/'n(fc±) 



Sn{k±) 



(fk dj ■ 17„(k) 



(10) 



(11) 



Here k± labels the slice and Sn{k±) is the region of the 
slice in which band n is occupied. Recalling Eq. (4) and 
noting that cij is the unit vector normal to the slice, the 
application of Stokes' theorem to Eq. (11) yields 



(l>n{k±) 



Cn(k±) 



A„(k) • dl 



(12) 



where Cn{kj_) is the oriented curve bounding Sn{k±) on 
the slice and (j)n{k±) has the interpretation of a Berry 
phase. For later convenience we also define 



<^(fcx)=^<^„(A:^) 



(13) 



and similarly Cj = ^„ c„j etc. The calculation of the 
AHC has thus been reduced to a calculation that is re- 
stricted to the Fermi surface only, in the spirit of Eq. (23) 
of Haldane's Ref. 10. 

In general, the occupied or imoccupied region of band 
n in slice k±_ need not be simply connected, in which case 
the boundary C„(fc_L) is really the union of several loops. 
Moreover, loops encircling hole pockets should be taken 
in the negative direction of circulation. This is illustrated 
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FIG. 1: Sketch of intersections of the Fermi surface with a 
constant-fcx plane. Open, hashed, and cross-hashed regions 

correspond to fiUing of zero, one, and two bands, respectively. 
The four small Fermi loops belong to the first band, while the 
large central one belongs to the second. Arrows indicate sense 
of circulation for performing the Berry-phase integration. 

in Fig. 1, where the first band exhibits four hole pockets 
and the second band has one electron pocket, so that 
Ci is the union of four countcrcirculating loops and C2 
is a fifth loop of positive circulation. If higher bands 
are unoccupied, then ^(fcx) for this slice is just given 
by the sum of the Berry phases of these five loops. We 
shall assume for simplicity in the following that Cn{k±) 
is simply connected, but the generalization to composite 
loops is straightforward. 



B. The quantum of Hall conductivity 

We claimed earlier that if band n is fully occupied, 

K„ in Eq. (7) is quantized to a reciprocal lattice vec- 
tor. This can now be seen by noting that under those 
circumstances the integral in Eq. (11) runs over a two- 
dimensional BZ, which can be regarded as a closed two- 
dimensional manifold (two- torus), and for topological 
roasons"'^^ the integral of the Berry curvature over such 
a closed manifold must be an integer multiple of 27r (the 
Chern number). Then each c„j is an integer, and K„ in 
Eq. (9) must be a reciprocal lattice vector as claimed. If 
the system is an insulator, then K = J2n (summed 
over occupied bands) is also guaranteed to be a recip- 
rocal lattice vector, and if it is a nonzero one, the in- 
sulator would have a quantized Hall conductivity and 
could be regarded as a quantum Hall crystal (or "Chern 
insulator" ).^°'^^ No physical realization of such a system 
is known experimentally, but the search for one remains 
an interesting challenge. 

Let us consider again a slice for which band n is fully 
occupied but has a non-zero Chern number. If this slice is 
regarded as an open rectangle (or parallelogram) rather 
than a closed two-torus, and a continuous choice of gauge 
is made in its interior (i.e., A„(k) is free of singularities), 
then the boimdary Cn{k±) is the perimeter of this rect- 
angle and Eq. (12) will yield the same integer multiple 



of 2tt as Eq. (11). In the spirit of Fig. 1, however, we 
prefer to regard the slice as a closed two-torus and to ex- 
clude the perimeter from our definition of the boundary 
Cn{k±). Then Cn{k±) is null and Eq. (12) vanishes for 
the case at hand, in disagreement with Eq. (11). The dis- 
agreement arises because of the impossibility of making 
a continuous choice of gauge on a closed manifold having 
a non-zero Chern number;^^ the best that can be done is 
to make A„(k) finite everywhere except at singularities 
("vortices") which, when included, restore the missing 
contributions of 27r. 

Returning to the general case of a partially occupied 
band n with Cn(k±) defined to exclude the perimeter of 
the slice, we conclude that Eq. (12) is really only guar- 
anteed to equal the true result of Eq. (11) modulo 27r. 
Moreover, the Berry phase will be evaluated in practice 
using a discretized Berry-phase formula^^ of the form 

<l>n{k±) = -Imln ]J(u„kj \unkj+i) (14) 
j 

where \cj discretizes the loop Cn{k±). (We will actu- 
ally use a modified version, Eq. (25), of this formula.) 
The choice of branch cut is now arbitrary, and again the 
agreement with Eq. (12) or Eq. (11) is only guaranteed 
modulo 27r. By convention one normally restricts phases 
to lie in the interval (— tt, tt], but then 4>n{k\_) would in 
general have unwanted clisc-ontinuities at some values of 
fc_L. In practice we discretize the k±_ integration, so that 
using Eq. (10), Cj = ^„ Cnj becomes 

1=1 

We then enforce continuity of the total phase <j){k±^ of 
Eq. (13) by choosing (^(i) such that \<j){i) — (f){i — l) \ <C 27r 
for each slice z = 2, 3, ... in sequence. Since the true phase 
given by the sum of contributions in Eq. (11) is also con- 
tinuous, this guarantees that our calculated (j){k±) dif- 
fers from the true one by the same multiple of 27r for 
all k±. Our computed AHC would then differ from the 
true one by a multiple of the quantum and could be said 
to give the "non- quantized part" of the intrinsic AHC in 
the sense of Haldane.^" However, it is straightforward to 
remove this overall ambiguity of branch choice by evalu- 
ating 4'{kj_) from Eq. (11) on the first slice and then en- 
forcing continuity for each subsequent slice, thus arriving 
at the correct AHC without any question of a quantum. 

We note in passing that an isolated point of degen- 
eracy ("Dirac point") between a pair of bands n and 
n+1 can generically occur in three-dimensional fc-space in 
the absence of time-reversal symmetry. If such a Dirac 
point occurs below the Fermi energy, then (f)n{k±) and 
(}>n+i{ki_) will, when evaluated from Eq. (11), exhibit 
equal and opposite discontinuities of 27r at the k± of the 
Dirac point. However, the total phase (f>{k±) will remain 
continuous, so that the algorithm described in the previ- 
ous paragraph will still work correctly. 
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We close this subsection by emphasizing that the dis- 
cussion of possible non-zero Chern numbers or the pres- 
ence of Dirac points is rather academic. In our calcula- 
tions on Fe, Ni and Co, we have not encountered any indi- 
cations of such anomalies; they presumably occur rarely 
or not at all in the materials studied here. 



C. Evaluation of the Fermi-loop Berry phase 

The essential problem now becomes the computation 
of the loop integral of Eq. (12). As is well known, the 
Berry connection A„(k) of Eq. (5) is gauge-dependent, 
i.e., sensitive to the fc-dependent choice of phase of the 
Bloch functions. If Eq. (12) is to be calculated by 
the direct evaluation of A„(k) and its subsequent inte- 
gration around the loop, this lack of gauge-invariance 
may present difficulties. For example, it means that 
there is no unique Kubo- formula expression for A„(k). 
An alternative and more promising approach is to com- 
pute (f)n{k±) by the discretized Berry-phase formula-' of 
Eq. (14), where the inner products are computed from 
the full first-principles calculations at neighboring pairs 
of fc-points around the loop. However, this may still be 
quite time-consuming if it has to be done at very many 
fc-points. We avoid this by using the technique of Wan- 
nier interpolation*'^' to perform the needed loop inte- 
gral inexpensively. In this formulation, the loop integral 
of Eq. (12) can be expressed as a sum of two terms, one 
in which a contribution to A„(k) is evaluated and inte- 
grated explicitly, and a second that takes a form like that 
of Eq. (14). 

The key idea of Wannier interpolation is to map the 
low-energy first-principles electronic structure onto an 
"exact tight-binding model" using a basis of appropri- 
ately constructed crystalline Wannier functions. For 
metallic systems like those considered here, the bands 
generated by these Wannier functions are only partially 
occupied. They are guaranteed by construction to repro- 
duce the true first-principles bands in an energy window 
extending somewhat above the Fermi level, so that all va- 
lence and Fermi-surface states are properly described.® In 
the Wannier representation, the desired quantities such 
as band energies, eigenstates and the derivatives of eigen- 
states with respect to wavevector k can then be eval- 
uated at arbitrary fc-points at very low computational 
cost. All that is needed is to evaluate, once and for all, 
the Wannier-basis matrix elements of the Hamiltonian 
and position operators.* It is worth pointing out that it 
may sometimes be expedient to drop some lower occupied 
bands and construct the Wannier functions so that they 
correctly represent the Bloch functions only in some nar- 
rower energy window containing the Fermi energy; since 
the present formulation involves only Fermi-surface prop- 
erties, the nonquantized part of the AHC will then still 
be given correctly. 

The Wannier construction procedure of Ref . 9 provides 
us with a set of M Wannier functions |Rn) (n = 1, M) 



in each cell labeled by lattice vector R. From these the 

Bloch basis functions \u'"^^) are constructed according to 
the Fourier transform relation 

\u^I^)=J2e-'^-^'^-^^\Rn). (16) 

R 

Here the superscript (W) indicates that these are ob- 
tained from the Wannier representation, that is, they are 
not yet Hamiltonian eigenstates. To obtain those, we 
construct the M x M Hamiltonian matrix 

HZ\^) = {u^I^\m\nT^) (17) 

via 

HiZ^=^e^^-^{On\H\Km). (18) 

R 

At any given k this matrix can be diagonalized to yield 
an M X M unitary matrix J7„m(k), i.e., 

?7t(k)ijW(k)?7(k) = ff(")(k) (19) 

where H^^\\i) = S^^Smn are the energy eigenvalues and 

l«ik') = El"-T>C^r„n(k) (20) 
jn 

are the corresponding band states. By the construction 
procedure of Ref. 9, is identical to the true (and 
similarly for the eigenvectors u'^) for all occupied states 
and low-lying empty states. This is strictly true only for 
fc-points on the original ab-initio mesh. The power of 
this interpolation scheme lies in the fact that, by virtue 
of the spatial localization of the Wannier functions, the 
error remains extremely small even for points away from 
that grid.^^ 

The next step is to evaluate on a two-dimensional 
mesh of fc-points covering a single slice and then use a 
contour-finding algorithm to map out and discretize the 
Fermi loops therein. This part of our scheme will be 
described in more detail in Sec. IIIB. For now we can 
just assume that the output is a sequence of points kj 
{j = 0, . . . , J — 1) providing a fairly dense mapping of 
the contour. (As before, we assume for simplicity that 
the Fermi contour consists of a single loop; the extension 
to multiple loops is straightforward.) 

Next we need to obtain the Berry connection A„ (k) = 

«^ikVkl4k^) as in Eq. (5). Using Eq. (20), this be- 
comes 

A4k) = J2uUk)A^\k)UmnO^) 

Im 

+iEc/t„(k)VkJ7„„(k) (21) 

m 

where 

A(r(k)=^^i^|Vk|e.^) (22) 
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is computed in practice from the expression 

Ai^^(k)=^e*-^(On|f|Rm) (23) 

R 

in a manner similar to Eq. (18). Details coneerning the 
method of calculating Eqs. (18) and (23) can be found in 
Ref. 8. 

The decomposition of A„(k) into two terms in Eq. (21) 
is an artifact of the choice of Wannier functions; only the 
sum of the two terms is physically meaningful (upon a cir- 
cuit integration). However, for a given choice of Wannier 

functions, the first term arises because the Bloch func- 

(a) 

tions !«,„!(. ) acquire some of the Berry curvature attached 
to the full subspace of M Wannier functions used to 
represent them, whereas the second term represents the 
Berry curvature arising from changes of character of this 
Bloch state within the Wannier subspace. To clarify this 
viewpoint, we introduce a notation^ in which ||u„k)) is 
defined to be the nth column vector of matrix U, so that 
the second term of Eq. (21) becomes i((i'nk|| VkH^^nk))- 
Plugging into Eq. (12), this yields 

Mi) = y((i;„k||AW(k)||t;„k))-dl 

+i ^{{vnk\\Vk\\vnu)) -dl (24) 

for the Berry phase of slice i appearing in Eq. (15). Note 
that the integrand in the first term is gauge- invariant 
(here "gauge" refers to the application of a phase twist 
ll^'rtk)) e^^^^^ ll^'nk))), while in the second term only 
the entire loop integral is gauge- invariant. Indeed, the 
second term is just a Berry phase defined within the M- 
dimensional "tight-binding space" provided by the Wan- 
nier functions. Recalling that kj for j = 0, . . . , J — 1 is 
our discretized description of the Fermi loop, and using 
standard methods for discretizing Berry phases^^ as in 
Eq. (14), our final result becomes 

j-i 

Mi) = ^((t;„k||A(w)(k)||t;„k)).Ak 
i=o 

,7-1 

-Imln []((t;„k,||t^„k,+i)), (25) 

where Ak = (kj+i — kj^i)/2. 

As we shall see below, in practice we only encounter 
closed orbits, in which case it is clearly appropriate to set 
kj = kg and close the phases with ||w„,k,,)) = H^'n.ko))- 
For lower-symmetry situations, however, open orbits 
with kj = ko -I- G may be encountered. Even in this 
case, however, we would still set Hiin^kj)) = Hfn.ko)); in 
contrast to the full Bloch states which obey^^ UnMj ~ 
e"'^'"" u„,ko) no extra phase factors are needed here be- 
cause the Fourier-transform convention of Eq. (16) treats 
the Wannier functions as though they are all nominally 
located at the cell origin. 



TABLE I: Calculated spin magnetic moment per atom (in 
/is) for the three transition metals Fe, Ni and Co, with mag- 
netization along [001], [111] and [001], respectively. 





bee Fe 


fee Ni 


hep Co 


Theory 


2.22 


0.62 


1.60 


Experiment"' 


2.13 


0.56 


1.59 



'Ref. 21. 



In summary, our strategy is to evaluate Eq. (15) by 
decomposing each generalized path C„(i) into connected 
simple loops, and sum the loop integrals as computed us- 
ing Eq. (25). The operations needed to evaluate Eq. (25) 
are inexpensive as they all involve vectors and matrices 
defined in the low- dimensional space of the Wannier rep- 
resentation. 



III. COMPUTATIONAL DETAILS 
A. First-principles calculations 

Fully relativistic band-structure calculations are car- 
ried out for the ferromagnetic transition metals Fe, Co 
and Ni at their experimental lattice constants (5.42, 
4.73, and 6.65 bohr, respectively) using the PWSCF 
codc.^*' Norm-conserving pseudopotentials with spin- 
orbit coupling^^ are generated using similar parameters 
as in Ref. 8. An energy cutoff of 60 Hartree is used 
for the planewave expansion of the valence wavefunc- 
tions (400 Hartree for the charge densities), and the 
PBE generalized-gradient approximation^* is used for 
the exchange-correlation functional. The self-consistent 
ground state is obtained using a 16 x 16 x 16 Monkhorst- 
Pack"'^^ mesh of fc-points and a fictitious Fermi smearing^'' 
of 0.02 Ry for the Brillouin-zone integration. 

The calculated spin magnetic moments are shown in 
Table I. The effect of spin-orbit coupling on these mo- 
ments is included in the calculation, since it is needed in 
any case to obtain a nonzero AHC. The agreement with 
experiment is rather good, confirming that our norm- 
conserving pseudopotentials are suitable for describing 
the ferromagnetic state of the transition metals. 

The maximally-localized Wannier functions are gen- 
erated using the WANNIER90 code;^^ details are given in 
Sees. IVA-IVC below. 



B. Mapping and sampling of Fermi loops 

As discussed above, our basic strategy involves divid- 
ing the BZ into a series of parallel slices and finding the 
intersections of the Fermi surface with each of these slices. 
Each slice is sampled on a uniform N x N fc-point mesh, 
with TV ranging from 300 to 500, and the band energies 
are computed on the mesh using Wannier interpolation. 
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FIG. 2: (Color online) Calculated Fermi-surface intersections 
(Fermi loops) on the kz = 0.02 plane for bee Fe; different 
bands are color-coded for clarity, (a) Fermi contours within 
the first Brillouin zone, (b) Fermi contours after reassembly 
to form closed contours by translating some portions by a 
reciprocal lattice vector. Inset: enlargement showing part of 
an avoided crossing where a refined mesh (black lines) is used 
to obtain a more accurate representation of the Fermi loop. 
The actual calculation is performed within the dashed box. 



A standard contour-finding algorithm of the kind used to 
make contour plots is then used to generate a hst of Fermi 
loops and, for each loop, a hst ko, . . . ,kj_i of fc-points 
providing a discretized representation of the loop. 

As shown in Fig. 2(a), the Fermi contours in the first 
BZ are sometimes composed of multiple segments termi- 
nating at the BZ boundary. To insure that we get closed 
loops suitable for the evaluation of Eq. (25), we actu- 
ally do the initial contour-finding procedure in an ex- 
tended zone with 3x3 times the size of the first BZ. We 
then select closed loops located near the central cell while 
identifying and discarding loops or portions of loops that 
correspond to periodic images of these chosen loops. The 
result is a set of closed loops that partially extend outside 
the first BZ as shown in Fig. 2(b). Of course, if there were 
open orbits on the Fermi surface, it would not always be 
possible to select closed loops in the above sense; one 
would have to accept a "loop" with kj = ko -I- G as dis- 
cussed following Eq. (25). However, we never encounter 



such open orbits in practice for the types of materials 
studied here, in which the magnetization is aligned with 
a three-fold, four-fold, or six-fold rotational symmetry 
axis. The slices are perpendicular to the symmetry axis, 
and the symmetry ensures that open orbits cannot occur 
on the slices. 

A potential difficulty in applying the Fermi-loop 
method to real materials arises from the possible pres- 
ence of degeneracies or near-degeneracies between bands. 
If two bands are degenerate at the Fermi energy, this 
means that two Fermi loops touch, and it is no longer 
straightforward to define and compute the Berry phases 
of these loops. Fortunately, the presence of ferromagnetic 
spin splitting and spin-orbit coupling removes almost all 
degeneracies. In our calculations we found no true de- 
generacies in hep Co or fee Ni, and the only degeneracies 
in bee Fe were found to lie in the = plane. (In the 
latter case, we avoid the = plane by picking a k± 
mesh that is offset so that this plane is skipped over.) 
On the other hand, we do find numerous weakly avoided 
crossings induced by the spin-orbit interaction, and while 
these introduce no difficulty in principle, they do require 
special care in practice. Indeed, we find that it is im- 
portant to sample the Fermi surface very accurately in 
the vicinity of these crossings. To do so, we calculate 
the Berry curvature at each kj using Wannier interpo- 
lation, and if a large value is encountered, we introduce 
a refined mesh with 4x4 greater density in this region, 
repeat the contour-finding procedure there, and replace 
the discretized representation of this portion of the loop 
with a denser one. We also take care to recompute SnW at 
each kj and iteratively adjust the fc-point location in the 
direction transverse to the loop in order to insure that 
£„k lies precisely at the Fermi energy. An example of 
a portion of a Fermi loop that has been refined in this 
way is illustrated in the inset to Fig. 2(b). Overall, the 
resulting number J of fc-points per loop ranges from sev- 
eral hundreds to thousands, depending on the size and 
complexity of the Fermi loop. 

In our current implementation, the entire procedure 
above is repeated independently on each of the slices. As 
already mentioned in Sec. II B, it is important to make 
a consistent choice of branch of the Berry phase on 
consecutive slices. We do this by adding or subtracting a 
multiple of 27r to the Berry phase calculated from Eq. (15) 
such that \(j)n{i) — — 1)1 ^ 27r is satisfied, always 
checking for consistency between the first and last slice. 



C. Use of symmetry to reduce computational load 

The presence of a net magnetization results in a con- 
siderable reduction in symmetry, but several symmetries 
still remain that can be exploited to reduce the computa- 
tional cost. In the previous Fermi-sea-based methods^'* 
the use of symmetries is straightforwardly implemented 
by restricting the fc-point sampling to the irreducible 
wedge of the BZ. For the Fermi-loop method, the use 
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FIG. 3: Illustration of use of and My mirror syrnrnctrios 
on a slice of the Brillouin zone in bcc Fe. Only the segment 
of the Fermi loop from Point 2 to Point 1 is actually com- 
puted; the three other segments are included using symmetry 
operations. 



of symmetries needs more careful treatment. 

Here we discuss the difficulties, and point out their 
solution, using ferromagnetic bcc Fe as an example. We 
focus our attention on the mirror symmetries and 

My. Since each slice lies in an x-y plane, we can use these 
to restrict the bandstructure calculation and the search 
for Fermi contours to a reduced BZ having one-fourth of 
the area of the full BZ, as shown by the dashed line in 
Fig. 2(b). However, a typical Fermi loop will no longer 
close within this reduced BZ. Because a Berry phase is a 
global property of a closed loop, one cannot just compute 
the Berry phase of open segment lying inside the reduced 
BZ and multiply by four; the Berry phase of this segment 
is ill-defined unless the phases of the wavefunctions at its 
terminal points are specified. 

Our solution to this difhculty is illustrated in Fig. 3. 
We make some arbitrary but definite choice of the phases 
of the Bloch functions in the upper-right segment, com- 
pute the open-path Berry phase following Eq. (25), and 
multiply by four. We then add corrections that take ac- 
count of the phase jumps at the segment boundaries. For 
example, we let Al^ acting on the Bloch states from 1 to 2 
define the Bloch states from 1' to 2'. The correction aris- 
ing from the I'-l boundary is then given by the phase of 
{uv\ui) = {MxUi\ui). (Here is defined in the spinor 
context and includes a complex conjugation component. 
Since the Bloch functions are expressed in the Wannier 
basis in our approach, information about the symmetries 
of the Wannier functions has to be extracted and made 
available for the application of the symmetry transforma- 
tions.) Similar corrections, using also My, are obtained 
for the 2'-2", 1"-!'", and 2"'-2 segment boundaries. By 
including these mismatch corrections, we are able to cal- 
culate the global Fermi-loop Berry phase in a correct and 
globally gauge-invariant manner. 

We have tested this procedure and confirmed that the 
results obtained are essentially identical to those com- 




FIG. 4: Calculated Fermi surfaces for bands 7-10 of bcc Fe 
(in order of upper left, upper right, lower left, lower right). 
The outside frame is the boundary of the Brillouin zone. 



puted without the use of symmetry. The BZ could in 

principle be reduced further in bcc Fe using the diagonal 
mirror operations, but we have not tried to implement 
this. 



IV. RESULTS 

In this section we present the results of our calculations 
of the anomalous Hall conductivity using the Fermi-loop 
approach of Eq. (25) as applied to the three ferromagnetic 
transition metals Fe, Co and Ni. 



A. bcc Fe 

We have previously presented calculations of the ARC 
of bcc Fe computed using the Fermi-sea formulation.^ 
Here we adopt the same choice of Wannier functions as 
in that work, namely 18 Wannier functions covering the 
s, p and d character and both spins. The orbitals of s, p, 
and Cg character are actually rehybridized into Wannier 
functions of sp^(P type, and the Wannier functions are 
only approximate spin eigenstates because of the pres- 
ence of spin-orbit interaction (see Ref. 8 for details). 

In our calculation for bee Fe, bands 5-10 cross the 
Fermi energy. Fig. 4 shows the Fermi-surface sheets 
for bands 7-10, plotted using the Xcrysden package. 
(Bands 5 and 6 give rise to small hole pockets, not 
shown.) Some of these sheets (especially 7 and 8) are 
quite complicated but, as expected, they all conform 
to the lattice symmetries. What is not clearly visible 
in these plots are the tiny spin-orbit-induced splittings, 
which change the connectivity of the Fermi surface. As 
mentioned earlier, such features play an important role 
in the AHC, and need to be treated with care. 
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FIG. 5: Calculated Berry phase <f){k±) of bcc Fe (in radians) 
as a function of k± (in units of 2-K/a). Solid line shows re- 
sults obtained from the Fermi-loop method of Eq. (25); circles 
indicate reference results obtained by the integration of the 
Berry curvature on each slice using Eq. (11). 



We take the magnetization to lie along the [001] axis. 
Choosing bi = (27r/o)(lT0) and bz = (27r/a)(110) in 
the notation of Sec. II A, it follows that as = 27rbi x 
b2/T4ccip — (0, 0, a) where Kecip is the primitive recipro- 
cal cell volume, and we only need to compute the c„3 in 
Eq. (10). The slices are square in shape, and k± = is 
discretized into 500 slices. 

In Fig. 5 we have plotted the total Berry phase Eq. (13) 
on each slice as computed from Eq. (25). The results are 
symmetric under mirror symmetry, so only half of the 
range of k±^ is shown. The sharp peaks and valleys in 
Fig. 5 are related to degenerate or near-degenerate bands 
that have been split by the spin-orbit interaction, as was 
illustrated, e.g., in the inset of Fig. 2. To validate the 
calculation, we compare it against a direct numerical in- 
tegration of the Berry curvature over the occupied bands 
using Eq. (11), as indicated by the symbols in Fig. 5. In 
spite of rather complex and irregular Fermi surfaces, the 
agreement between the two methods in Fig. 5 is excellent. 

The values of the integrated anomalous Hall conduc- 
tivity using the new approach and the reference approach 
are shown in the first and third lines of Table II. The sec- 
ond line shows the contribution obtained from integrating 
only the first term of Eq. (25); clearly, this contribution 
is very small. The agreement with the previous theory 
of Yao et alJ is excellent, while the agreement with ex- 
periment is only fair. Table II will be discussed further 
in Sec. IV D. 

Our approach also opens the possibility of discussing 
which Fermi sheets are responsible for features visible in 
Fig. 5. For example, the dip near k± = 0.03 and the 
peak near k± =0.33 (in units of 27r/a) come from sheets 
8-9, the peak near k± = 0.09 comes from sheets 6-8, and 
the complex structure in the range of k± from 0.36-0.50 
comes mainly from sheets 7-9. Overall, the contribution 
from bands 5 and 10 are almost negligible, and bands 7-9 
give much the largest contributions. 



TABLE II: Anomalous Hall conductivity, in S/cm. First three 
rows show values computed using Eqs. (6)-(10) together with 
Eq. (25), the first term only of Eq. (25), or Eq. (11), respec- 
tively. Results of previous theory and experiment are included 
for comparison. 





bcc Fe 


fee Ni 


hep Co 


Fermi loop 


750 


-2275 


478 


Fermi loop (1st term) 


7 





-4 


Berry curvature 


753 


-2203 


477 


Previous theory 


751" 


-2073*" 


492'' 


Experiment 


1032= 


-646"^ 


480*^ 



''Ref. 7. 

''Y. Yao, private communication. 
=Ref. 24. 
•^Ref. 25. 
=Ref. 26. 



B. fee Ni 

For fee Ni we chose 14 Wannier functions, seven each of 
approximately spin-up and spin-down character. These 
were comprised of five Wannier functions of d-like sym- 
metry centered on the Ni atoms and two Wannier func- 
tions of tetrahedral symmetry located on the tetrahedral 
interstitial sites, similar to the choice that was made for 
Cu in Ref. 9. The inner energy window was chosen to 
extend 21 eV above the bottom of the bands, thus ex- 
tending 7.1 eV above the Fermi energy and including 
several unoccupied bands as well. 

Our calculation is consistent with previous DFT cal- 
culations in predicting that five bands (bands 8-12) cross 
the Fermi energy in fee Ni. The Fermi sheets for bands 
9-12 are shown in Fig. 6. Band 8 only barely crosses the 
Fermi energy and gives rise to very small hole pockets 
near the X points (even smaller than those illustrated for 
band 9). The existence of these pockets is a delicate fea- 
ture that is not clearly confirmed experimentally and is 
inconsistent with some recent LDA + U calculations.^^ 
However, including them or not has very little influence 
on our calculated ARC, as explained below. The shapes 
of the Fermi sheets in fee Ni are somewhat more spherical 
than those of bcc Fe. As expected, they again conform 
to the lattice symmetries. 

In the case of fee Ni, the magnetization lies along 
the [111] axis. Choosing bi = (27r/a)(022) and b2 = 
(27r/a)(202) in the notation of See. II A, it follows that 
as = 27rbi x b2/K-ecip = {a, a, a) = a"\/3e(iii), and we 
only need to compute the c„3 in Eq. (10). The slices are 
hexagonal in shape, and k± = k -6(111) is discretized into 
about 100 slices. 

The results are plotted in Fig. 7, along with symbols 
denoting the reference calculation by an integration of the 
Berry curvature over the slice. Once again, the agree- 
ment is very satisfactory. The values of the integrated 
AHC are again summarized in Table II. A band-by-band 
analysis indicates that band 8 gives only a very small 
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FIG. 6: Calculated Fermi surfaces for bands 9-12 of fee Ni. 
The outside frame is the boundary of the Brillouin zone. 



FIG. 8: Calculated Fermi surfaces for bands 18-21 of hep Co. 
The outside frame is the boundary of the Brillouin zone. 



contribution, less than 5% in magnitude and opposite in 
sign, to the total AHC. The hole pockets in band 9 give a 
slightly larger positive contribution, but we find that the 
dominant negative contribution to the AHC comes from 
bands 10-12. 



C. hep Co 

Co in the hep structure has two atoms per unit cell. 

We choose 18 Wannier functions per Co atom, nine for 
each spin, in a very similar manner as was done for Fe in 
Sec. IV A. We therefore have 36 Wannier functions per 
cell. 
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FIG. 7: Calculated Berry phase <j}{k±) of fee Ni (in radians) 
as a function of k± (in units of 2-K/\/3a). Solid line shows 
results obtained from the Fermi- loop method of Eq. (25); cir- 
cles indicate reference results obtained by the integration of 
the Berry curvature on each slice using Eq. (11). 



In our calculation, seven bands (bands 16-22) cross the 
Fermi energy in hep Co. We show the four largest Fermi- 
surface sheets associated with bands 18-21 in Fig. 8. The 
Fermi surfaces can be seen to respect the 6-fold crystal 
symmetry, and none of them touch each other. 

The magnetization of hep Co lies along the [001] axis. 
We thus choose bi = (27r/a)(l/-\/3,— 1,0) and b2 = 
(27r/a)(l/-\/3, 1,0) in the notation of Sec. II A, and it fol- 
lows that as = 27rbi X b2/14ecip = (0,0, c). The slices 
are hexagonal in shape, and k± = is discretized into 
about 200 slices. 

The results are plotted in Fig. 9, along with the sym- 
bols denoting the reference calculation by integration of 
the Berry curvature. Once again, the peaks and valleys 
correspond to the places where two loops approach one 
another closely. Some pieces of the Fermi surfaces of hep 
Co are nearly parallel to the slices (see the bottom right 
panel of Fig. 8), so that the number and shapes of the 
Fermi loops sometimes change rapidly from one slice to 
another. In particular, we found it difficult to enforce 
continuity of the branch choice of Eq. (25) as a func- 
tion of k±_ near the sharp features at ki_a/2-n = 0.18 and 
0.42 in Fig. 9. We therefore redetermined the correct 
branch choice by comparing with the result of the Berry- 
curvature integration at slices just outside these difficult 
regions. Despite these difficulties, it can still be seen that 
the Fermi-loop method works well for this case. Some of 
the sharp structure appearing in Fig. 9 in the range of fcj, 
from 0.4-0.5 (in units of 2ti / a) arises from the small hole 
pocket in band 16, but this gives a rather small contribu- 
tion to the total AHC. The peak around k±_ = 0.14 and 
the sharp dip around 0.18 comes mainly from the sheets 
associated with bands 20 and 21. 
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FIG. 9: Calculated Berry phase (l>{k±) of hep Co (in radians) 

as a ftmction of k± (in units of 2n/c). Solid line shows results 
obtained from the Fernii-loop method of Eq. (25); cireles indi- 
cate referenee results obtained by the integration of the Berry 
eurvature on each slice using Eq. (11). 



D. Discussion 

1. Internal consistency of the theory 

The second row of Table II shows the results computed 
using only the first term of Eq. (25) . In each case, its con- 
tribution is less than 1% of the total, and would therefore 
be negligible for most purposes. Actually, it can be shown 
that the inclusion of the first term only in Eq. (25) of the 
present method is equivalent to carrying out the Berry- 
curvature integration approach of Ref. 8 with the D D 
term omitted in Eq. (32) of that work (that is, only the 
D-A and Q terms included). We have carried out this 
comparison and find values of 7, —0.5 and — 2S/cm for 
bcc Fe, fee Ni, and hep Co, respectively, in very good 
agreement with the values reported in Table II. The 
physical interpretation for the small terms in the second 
row of Table II is basically that the full set of Bloch-like 
states constructed from the Wannier functions (e.g., the 
manifold of 18 Bloch-like states in bcc Fe) has some small 
Berry curvature of its own, and the projection of this cur- 
vature onto the occupied subspacc gives the small first 
term of Eq. (25). On the other hand, spin-orbit induced 
splittings across the Fermi level between Bloch-like states 
built from these Wannier functions give large, sharply 
peaked contributions to the Berry curvature of the occu- 
pied subspace, and make a very much larger contribution 
to the total AHC. Of course, the precise decomposition 
between the first and second term of Eq. (25) depends on 
the exact choice of Wannier functions, but the present re- 
sults seem to indicate that the dominance of the second 
term is probably a general feature, at least for systems in 
which the Wannier functions are well localized and the 
spin-orbit splitting is not very strong. 

As mentioned in the previous section, the overall agree- 
ment seen in Table II between the results computed us- 



ing the Fcrmi-loop approach and those computed using 
the Berry-curvature integration indicate the internal con- 
sistency of our theory and implementation. The agree- 
ment with the results of Yao and coworkers, which were 
obtained by a Berry-curvature integration using an all- 
electron approach,'' also demonstrates the robustness of 
our pseudopotential implementation, including its ability 
to represent spin-orbit interactions correctly. 



2. Comparison with experiment 

In the last row of Table II we show comparison with 
some representative experimental values for the AHC of 
Fe, Ni, and Co. However, it should be kept in mind that 
there is some uncertainty and variation in the values re- 
ported by different groups. For example, Ref. 28 gives a 
value for Ni of — 753 S/cm and Ref. 29 reports a value for 
Co of 500 S/cm. It could well be that different kinds of 
experimental samples have different impurity and defect 
populations, leading to different extrinsic contributions 
to the AHC. Since the theoretical values presented in 
Table II are all computed by including only the intrinsic 
Karplus-Luttinger contribution to the AHC, so that ex- 
trinsic skew and side-jump scattering contributions are 
neglected, it is most appropriate to compare with experi- 
mental measurements in which the effects of the intrinsic 
contribution are isolated. 

A serious effort in this direction has recently been made 
by studying and correlating the variation of both the 
longitudinal and the anomalous Hall conductivity as a 
function of temperature.^^ It was found that for Fe, Co, 
and Ni, \axy \ remains roughly constant between 150 and 
300 K while a^x changes by about a factor of four. The 
value of axy in this plateau was attributed to the intrinsic 
mechanism, which should be independent of the scatter- 
ing rate. The values thus obtained are about 970 and 
—480 S/cm for Fe and Ni films respectively (the value 
quoted in Table II for Co from the same work is also a 
film value), and about 2000 S/cm for single-crystal Fc. 
The factor-of-two difference reported in Ref. 26 between 
the intrinsic axy of Fe in the single-crystal and film forms 
is puzzling and deserves further investigation. 

Turning now to the comparison between theory and ex- 
periment, we find a very rough agreement at the level of 
signs and general trends. However, the agreement is not 
quantitatively accurate, except for Co where the agree- 
ment is good. For Fe our results arc in very rough agree- 
ment (~25% low) compared to the results of Ref. 24 or 
the film results of Ref. 26, but a factor of two smaller 
than the single-crystal results of Ref. 26. Clearly the 
most serious discrepancy is for Ni, for which we get a 
consistent sign but a much larger magnitude than indi- 
cated by the experiments. (The issue^"^ of whether a hole 
pocket really appears in band 8 of fee Ni is not relevant 
since, as indicated in Sec. IV B, it makes a quite small nu- 
merical contribution to our theoretical result.) While the 
available experimental values for Ni appear to be roughly 
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consistent with each other, we are not aware of any study 
using the methods of Ref. 26 apphed to both films and 
single crystals of Ni. Until the experimental situation is 
clarified further, a final judgment on the degree of dis- 
agreement with the values based on DFT calculations 
should perhaps be withheld. 

In the meantime, it would be desirable to face some of 
the challenges and open questions that remain on the the- 
oretical side. For example, not much is yet known about 
the accuracy of common exchange-correlation function- 
al, such as the PBE functional used here,^^ for com- 
puting the AHC. Fe, Ni and Co have open d shells and 
can be considered from one point of view to be strongly- 
correlated systems. The fact that the magnetic mo- 
ments are given accurately by DFT (sec Table I) does 
not necessarily mean that more delicate properties, es- 
pecially those like the AHC and the magnetocrystalline 
anisotropy^^ that depend on spin-orbit interactions, will 
be given accurately. It is possible that the use of more so- 
phisticated density- functional approaches (e.g., current- 
density functional theories) or higher-level many-body 
approaches may ultimately prove necessary. Finally, it 
would be desirable to develop DFT-based methods for 
computing the defect-related extrinsic contributions, but 
this will also prove to be a daunting challenge, not least 
because the relevant defect populations are not known. 

In summary, experiments and DFT-based theories 
agree on the orders of magnitude and signs of the intrin- 
sic Karplus-Luttinger contributions to the AHC in these 
three ferromagnetic metals, and the results for Fe and es- 
pecially Co suggest that quantitative agreement may be 
obtainable. The substantial discrepancy for Ni deserves 
further attention. 



V. COMPUTATIONAL EFFICIENCY 

The motivation for developing a method for computing 
the AHC that relies only on information computed on the 

Fermi surface is, to some degree, esthetic and philosoph- 
ical: Haldane argued that the AHC is physically most 
naturally regarded as a Fermi-surface property,^" and as 
such should be computed using a method that does not 
make use of extraneous information in arriving at the de- 
sired quantity. However, a much more important motiva- 
tion from the practical point of view is the idea that the 
computational effort might be drastically reduced by hav- 
ing to compute quantities only on the two-dimensional 
Fermi surface rather than on a three-dimensional mesh 
of fc-points. 

In the present implementation as it stands, unfortu- 
nately, the computational savings gained through the use 
of the Fermi-loop Berry-phase approach is quite modest. 
After taking advantage of the symmetry as discussed in 
Sec. Ill C, the total computational time of our AHC cal- 
culation for bcc Fe is about 1.7 hours using a 200 x 200 
k-mesh on each of 500 slices, to be compared with about 
2 hours using our previous method of Ref. 8. (These 



timings are on a 2.2 GHz AMD-Opteron PC, and nei- 
ther includes the Wannier construction step, which takes 
about 2.5 hours.) Roughly, the work on each slice can 
be divided into three phases: Step 1, computing the en- 
ergy eigenvalues on the 200 x 200 k-mesh; Step 2, execut- 
ing the contour-finding algorithm; and Step 3, evaluating 
Eq. (25) on the discretized Fermi loops. We find that less 
than 1% of the computer time goes to Step 2, while the 
remainder is roughly equally split between Step 1 and 3. 
The operations in these steps have been greatly acceler- 
ated by making use of Wannier interpolation methods, 
but this is also the case for the comparison method of 
Ref. 8. (We emphasize that, for this reason, both the 
method of Ref. 8 and the present one are orders of magni- 
tude faster than methods based on direct first-principles 
calculations at every /c-point.) 

Many opportunities for further reduction of the com- 
puter time are worthy of further exploration. Regarding 
Step 1, for example, at the moment the contour-finding is 
done independently on each slice; it might be much more 
efficient to step from slice to slice and use a local algo- 
rithm to determine the deformation of the Fermi contours 
on each step. It may also be possible to do a first cut at 
the contour-finding using a coarser k-mesh (say 50 x 50) 
and then refine it in regions where the loops approach 
one another or have sharp bends. It may also be possible 
to take larger steps between slices in most regions of k±_ , 
and fall back to fine slices only in delicate regions. In 
implementing all such strategies, however, one should be 
careful to avoid missing any small loops that might ap- 
pear suddenly from one slice to the next, or which might 
be missed on an initial coarse sampling of the slice. It 
may also be interesting to explore truly three-dimensional 
algorithms for finding contour surfaces, and then derive 
two-dimensional loops from these. 

As for Step 3, it should be possible to use a lower 
density of fc-points in the portions of the loop discretiza- 
tion where the character of the wavefunctions is changing 
slowly. The time for this step will also obviously benefit 
from taking larger steps between slices in regions where 
this is possible. Finally, a reduction by a factor of two 
or more may be possible by making use of symmetries 
not considered in Sec. HI C, such as the diagonal mirror 
symmetries {x y etc.) in bcc Fe. 

The exploration of these issues is somewhat indepen- 
dent from the quantum-mechanical formulation of the 
underlying theory, which is the main focus of the present 
work, and wc have therefore left the exploration of these 
possibilities for future investigations. 

Finally, it should be emphasized that the computa- 
tional load scales strongly with the dimension of the Wan- 
nier space used to represent the wavefunctions. In our 
calculations, this was 18, 14, and 36 for Fe, Ni, and Co, 
respectively. In some materials, there may be only a few 
bands crossing the Fermi energy, and it might be possi- 
ble to represent them using a much smaller number of 
Wannier functions. This is the case in many transition- 
metal oxides such as Sr2Ru04, cuprate superconductors. 



12 



etc. In ferromagnetic materials of this kind, it should 
be possible to choose an inner window in the Wannier 
disentanglement procedure^ that brackets the Fermi en- 
ergy but does not extend to the bottom of the occupied 
valence band, and to generate just a handful of Wannier 
functions (e.g., three t2g orbitals times two for spin) to 
be used in the Wannier interpolation procedure. Then 
all matrices used in that procedure would be very much 
smaller (e.g., 6x6) and the computation would go con- 
siderably faster. 

VI. SUMMARY 

In summary, we have developed a first-principles 
method for computing the intrinsic AHC of ferromagnets 
as a Fermi-surface property. Unlike conventional meth- 
ods that are based on a A:-space volume integration of the 
Berry curvature over the occupied Fermi sea, our method 
implements the Fermi-surface philosophy by dividing the 
Brillouin zone into slices normal to the magnetization 
direction and computing the Berry phases of the Fermi 
loops on these slices. While Haldane has pointed out 
that only the non-quantized part of the AHC can be de- 
termined in principle from a knowledge of Fermi-surface 
properties only, we find in practice that it is straight- 
forward to make the correct branch choice and resolve 



the quantum of uncertainty by doing a two-dimensional 
Berry-curvature integration on just one or a few of the 
slices. Our method also makes use of methods of Wannier 
interpolation to minimize the number of calculations that 
have to be done using a full first-principles implementa- 
tion; almost all the operations needed to compute the 
AHC are actually done by working with small matrices 
(e.g., 18x18 for bee Fe) in the Wannier representation. 
The new method also allows us to disciiss the contribu- 
tions to the AHC arising from individual Fermi sheets or 
groups of sheets. 

We have tested and validated our new method by 
comparing with our earlier implementation of a Fermi- 
sea Berry-curvature integration for bcc Fe, fee Ni, and 
hep Co. The different crystal structures and magnetiza- 
tion orientations in these three materials also allow us 
to demonstrate the flexibility of the method in dealing 
with these different cases. We find excellent agreement 
between the two approaches in all cases. 
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